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ABSTRACT: The /iz/SSM provides a solution to the \x problem of the MSSM and explains 
the origin of neutrino masses by simply using right-handed neutrino superfields. We have 
completed the analysis of the vacua in this model, studying the possibility of spontaneous 
CP violation through complex Higgs and sneutrino vacuum expectation values. As a conse- 
quence of this process, a complex MNS matrix can be present. Besides, we have discussed 
the neutrino physics and the associated electroweak seesaw mechanism in the //z/SSM, in- 
cluding also phases. Current data on neutrino masses and mixing angles can easily be 
reproduced. 
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1. Introduction 

Although the minimal supersymmetric extension of the standard model (MSSM) 
reveals as a solution to the hierarchy problem, we still remain puzzled about the 
origin of the /i-term in the superpotential, known as the /^-problem |l[]. On the 
other hand, the fact that neutrinos are not massless [0 suggests that the MSSM is 
incomplete. Motivated by these two facts, the " [i from v n supersymmetric standard 
model (/iz/SSM) |3], H], which relies on the existence of right-handed neutrinos, 
arises as an alternative to the MSSM, providing a solution to the /i-problem and 
explaining the origin of neutrino masses. 

In particular, the superpotential of the /iz/SSM contains, in addition to the usual 
Yukawa couplings for quarks and charged leptons, Yukawa couplings for neutrinos 
H u L z> c , terms of the type £> c HdH u producing an effective fi term through right- 
handed sneutrino vacuum expectation values (VEVs), and also terms of the type 
i> c v c i> c avoiding the existence of a Goldstone boson and contributing to generate 
effective Majorana masses for neutrinos at the electroweak scale. Actually, the ex- 
plicit breaking of R-parity in this model by the above terms produces the mixing of 
neutralinos with left- and right-handed neutrinos, and as a consequence a general- 
ized matrix of the seesaw type that gives rise at tree level to three light eigenvalues 
corresponding to neutrino masses Hj. 
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Following this proposal, several papers have studied different aspects of the 
yUz/SSM. In H the parameter space of the /ii/SSM was analyzed in detail, study- 
ing the viable regions which avoid false minima and tachyons, as well as fulfill the 
Landau pole constraint. The structure of the mass matrices, and the associated 
particle spectrum was also computed, paying special attention to the mass of the 
lightest Higgs. In || neutrino masses and mixing angles were discussed, as well as 
the decays of the lightest neutralino to two body (WMepton) final states. The corre- 
lations of the decay branching ratios with the neutrino mixing angles were studied as 
another possible test of the /i^SSM at the LHC. The phenomenology of the /iz/SSM 
was also studied in @, particularized for one and two generations of right-handed 
sneutrinos, and taking into account all possible final states when studying the decays 
of the lightest neutralino. Possible signatures that might allow to distinguish this 
model from other R-parity breaking models were discussed qualitatively in the last 
two papers. Let us finally mention that terms of the type v c HdH u and i> c v c v c were 
also analysed as sources of the observed baryon asymmetry in the Universe || and 
of neutrino masses and bilarge mixing ||, respectively. 

The goal of this work is twofold; first, we complete the analysis of the vacua of 
the /zz/SSM presented in 0, studying spontaneous CP violation (SCPV) of the tree- 
level neutral scalar potential. In particular, we explore CP violation in the lepton 
sector and show how phases for the tree-level Maki-Nakagawa-Sakata matrix (MNS) 



llj may arise due to the fact that the minimum of the scalar potential with real 
parameters has complex VEV solutions. Second, we discuss neutrino physics and the 
seesaw mechanism in the /iz/SSM, including also phases. 

Let us recall that, although there is evidence for CP violation in the quark sector 
of the standard model, there are not experimental traces of it in the leptonic part. CP 
can be explicitly broken through complex parameters in the Lagrangian or can arise 
spontaneously in a CP conserving Lagrangian (e.g. with all the parameters being 
real) through complex VEVs. Although the standard model as well as the MSSM 
do not allow for SCPV, in more complicated models both sources of CP violation, 
complex parameters and complex VEVs, could be present. 

Concerning the quark sector, a recent study argues that the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix is likely complex |TT|. This conclusion is supported by the 
measurement of the unitarity triangle angle 7 by BaBar and Belle collaborations 
JT2] , P~5| . This evidence of a complex CKM matrix has ruled out Next-to-MSSM 
(NMSSM)-like models with SCPV (see e.g. [£L4[]) for being the entire source of CP 
violation in the quark sector, since the CKM matrix in such models is real. Thus 
complex parameters are necessary in the quark sector. Given the structure of the 
/iz/SSM, this fact also holds for this model. On the other hand, as mentioned above, 
we will show that SCPV can be generated in the leptonic sector of the /xz/SSM, as 
well as phases for the MNS matrix. 

One argument in favor of the presence of SCPV at the Lagrangian level is that, 
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if the determinant of the quark mass matrix is real, it leads to a solution to the 



strong CP problem [15|. Extensions of the MSSM having this property, have been 



extensively studied in the literature (see e.g. In those scenarios, the quark 

sector of the model is extended in such a way that the effective 3x3 CKM matrix 
is complex whereas the determinant of the quark matrix is real. 

Other authors have extended the Higgs sector of the models, leading to SCPV 
with a complex CKM matrix |T7]]. Last but not least, in supersymmetric (SUSY) 
models with both CP and Peccei-Quinn symmetries, SCPV can be used as a solution 
to the SUSY phase problem |fl8[| . 

Regarding extensions of the /iz/SSM, the SCPV scenario with a complex CKM 
matrix can be accomplished by adding two more families of Higgs doublets. In this 
case the model would contain three families of matter and Higgs fields. This pos- 
sibility is well motivated phenomenologically, since the potential problem of flavor 
changing neutral currents can be avoided |19|]. In addition, having three Higgs fam- 



ilies is favored in some string scenarios [20]. Indeed, extensions of the quark sector 



of the model can also be studied, without altering the results here presented. 

What we want to point out in this work is that SCPV is possible in the simplest 
version of the /i^SSM, i.e. with only one family of Higgs doublets, and therefore it is 
worth studying its consequences. Following this philosophy, the paper is organized 
as follows. Section |2| is devoted to complete the analysis of the vacuum of the 
/iz/SSM started in including SCPV solutions. In Section ^|we examine the seesaw 
mechanism as the origin of neutrino masses and mixing angles in the model. In 
Section [| we carry out a detailed numerical analysis of the tree-level neutral scalar 
potential, showing explicitly that SCPV solutions are possible, and discussing their 
implications on the neutrino sector of the model. Finally, the conclusions are left 
for Section |5|. Minimization equations of the model and an approximate analytical 
formula for neutrino masses are given in the Appendices. 



2. Complex VEVs in the /iz/SSM 

The superpotential of the /zz/SSM introduced in || is given by 

w = E h *i & ^ + Y ^ ^ <K % + Y ^ ^ Lb i ^ + Y ^ *i l * »i) _ 

o,6 i,j 

-EE ^ % v a d H b u + J2 \^tm > (2- 1 ) 

a,b i i,j,k 

where we take Hj = (H d ,H~), E T U = (H+,H°), Qf = Lj = (i>,,e L J, 

i,j, k = 1,2,3 are family indices, the 3x3 matrices Y are dimensionless Yukawa 
couplings, a, b = 1, 2 are SU{2)l indices and ei2 = 1. As mentioned in the Introduc- 
tion, in addition to the MSSM Yukawa couplings for quarks and charged leptons, the 
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/iz/SSM superpotential contains Yukawa couplings for neutrinos, and two additional 
type of terms involving the Higgs doublet superfields, H d and H u and the three right- 
handed neutrino superfields, Of, with the dimensionless vector coupling A and the 
totally symmetric tensor k. 

As discussed in [[|, when the scalar components of the superfields z> c , denoted 
by £>f, acquire VEVs of the order of the electroweak scale, an effective interaction 
\xH\H2 is generated through the fifth term in Eq. (p.l|) , with fi = Aj(z>f). The last 
type of terms in Eq. (|2.1|) is allowed by all symmetries, and avoids the presence of an 
unacceptable Goldstone boson associated to a global U(l) symmetry. In addition, 
it generates effective Majorana masses for neutrinos at the electroweak scale. These 
two type of terms break explicitly .R-parity and lepton number. 

Working in the framework of gravity mediated supersymmetry breaking, the 
Lagrangian £ soft is given by: 



-C 



soft 



1,3 



c ij J ij J 



+ E K< m* m + ™k K*K] 

a 

+ E E tab [i A - Y -h H u Qi ^ + ( A dY d )v H a d Q\ d] + (A e Y e )ij H a d L 



1 c j 



a,b i,j 



+ (A 1/ Y 1/ ) ij H b u L?u° + c.c. 



+ 



E E ^{Ax\)i vf H a d H b u + \{A K K) ijk vtv'pi + c.c. 



ah i 



M 3 A 3 A 3 + M 2 A 2 A 2 + Mi Ai Ai + c.c. 



(2.2) 



In addition to terms from £ so ft, the tree-level scalar potential receives the D and 
F term contributions also computed in ||. In the following we will suppose that 
CP is a good symmetry of the model, taking all the parameters in the neutral scalar 
potential real and assuming that CP is only violated by the VEVs of the scalar fields 

(H° d ) = e^v d , (H° u ) = e^v u , (v i ) = e^v i , (pfi = e*f v\ . (2.3) 

We then obtain for the tree-level neutral scalar potential, 



V° = Koft + V D + V F 



(2.4) 
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where 



V so ft = m 2 Hd v d v d + m 2 Hu v u v u + Y]m\ mm v { Vj cos(x* - Xj) + Y] ™Jp.^, c cos(<p„c - tpj) 



+ 2 ^2(AvY„)ijV u iSiiSj cos(xi + y v <r) , 



(2.5) 



G 



(2.6) 



with G 2 = g 2 + gj, and 

W = £(A,) 2 « + 

+ ^ C0S(^ - + ^ Wvl^ COs(^c - y^c) 

i,j,k,l m 

+ 2 ~y^y^ K ikj \ k V d V u VtVj COs((/?„c + (f u c - 

+ ^ Y >> j i K iikV»i' j vfuZ cos(^c + y?^ - Xj) 
- 5^ 5^ ^Aj-Ud^i/i cos(</^ - Xi 

i j 

+ X! X! Y »ik Y » i k v l 1 'i l 'i cos (Xi - Xj) 

i,j k 

+ J] ^^^^^cos^ - ^c) . 



(2.7) 



We observe that in the potential there are seven independent phases, and we have 
defined them as 

<Pv = tPv u + <Pv d , Xi = + ¥v u , Vvf- (2.8) 

Now one can derive the fifteen minimization conditions with respect to the mod- 
uli v d , v u , z/ c , z/j, and phases <p v , Xh an d Vu c ■ These are written in Appendix A. 
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Finding minima requires the solutions of equations (|A.1HA.7|) . A standard way to 
obtain this is to give the values of the cosines of the phases in terms of the mod- 



uli, using the triangle method [21, 22, 23] for the equation of the phases, and then 



substitute the expressions in the minimum equations for the moduli, solving them 
numerically. This method permits to demonstrate the existence at tree level of only 
real minima in several models. This is for example the case of the NMSSM [ 23] , and 



the MSSM with extra doublets. The latter result has been proved for the MSSM with 
an extra pair of Higgs doublets |22| (the so called 4D model), the bilinear R-parity 
violation model (analogous to a 5D model because of the VEVs of the left-handed 
sneutrinos), and the MSSM with two extra pair of Higgs doublets (6D model) [p3| . 

Another way of finding minima consists of using as input the phases and solve 
the fifteen equations to fix the variables that are linear in these equations, as it is the 
case of some of the soft terms. This is the procedure that we will follow in Section 4. 

A simple way to prove the existence of CP violating minima in the /zi^SSM 
is using the results of Ref. p3| , where the authors prove that SUSY scenarios for 



SCPV require singlets. In particular, they found that, if the singlets do not introduce 
dimensional parameters in the superpotential (i.e. no linear or bilinear terms), the 
MSSM extended with two gauge singlets would be the minimal SUSY model where 
CP violation can be generated spontaneously. Since that model is a limiting case of 
the /xz/SSM with vanishing neutrino Yukawa couplings Y Vi . = 0, A3 = 0, and k 333 = 
^322 = ^332 = ^311 = ^331 = «i23 = 0, this would prove that the /zz/SSM can break 
CP spontaneously. Let us remark that, since in the /xz/SSM one is using a seesaw 
at the electroweak scale, the Y v .. have to be very small compared with the other 
parameters ||, and as a consequence the neutral scalar potential can be understood 
as a deformation of the MSSM extended with three gauge singlets. Although there is 
no literature about general solutions that break CP spontaneously in the latter, it is 
obvious that this model contains the MSSM extended with two singlets as a limiting 
case when k 333 = k 322 = k 33 2 = «3ii = ^331 = ^123 = 0, and A 3 = 0. As already 



mentioned, SCPV solutions are well known in this case |23], 25|. Thus one could 



argue that a subset of solutions with neutrino masses different from zero could be 
obtained deforming the scalar potential of the MSSM extended with three singlets 1 
through non-zero Y n .. 

In Sect. |]we will do a thorough numerical analysis showing explicitly how SCPV 
is realizable in the leptonic sector. Nevertheless, it is worth pointing out here that to 
find complex solutions is a non-trivial task compared to the search of real ones. As 
we will show, the key of SCPV is on the (A K K)ijk terms used as inputs. In order to 



1 Since only mass differences for neutrino masses have been measured, in principle two right- 
handed neutrino supermultiplets are enough to give two tree-level masses and also break CP spon- 
taneously. Thus a version of the /ij/SSM with only two right-handed neutrinos instead of three could 
be formulated. Nevertheless, we will follow the philosophy that the existence of three generations 
of all kind of leptons is more natural. 
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fulfill the minimization equations, the basic requirement is that entries different from 
(A K K,)m must be allowed. In addition, these parameters have to be chosen carefully 
to obtain SCPV as a global minimum. 

In the next section we will study the seesaw mechanism in the model as the 
origin of neutrino masses and mixing angles. 



3. Neutrino masses and mixing angles 

In the /zi/SSM the MSSM neutralinos mix with the left- and right-handed neutrinos as 
a consequence of R-parity violation. Therefore the right-handed neutrinos behave as 
singlino components of the neutralinos. In the basis x° T — (B°, W°, H d , H u , v Ri) v L .) 
the neutralino-neutrino mass matrix was given in [Q, |J for real VEVs. Considering 
now the possibility of complex VEVs the result is given by 



M, 



M m 
m T 3x3 



(3.1) 



where the neutralino mass matrix is 



M 



Mi 


-A(H° d y 

A(H°uT 










M 2 

B(H° d y 

-B(H°)* 






-A(H° d )* 
B(H° d )* 


-A3<flS> 



A(H°)* 
-B(H° U )* 

-Ai<£f) 


-\ x {H°) + Y Vi ,(n) 

-\ 2 (H° d ) + Y Ui2 (9 t ) 
-X 3 (H° d ) + Y Vi3 {Di) 






-Ai(fl2> 







MH° d ) + Y Vi2 (n 

2Kl2j(l>j} 
2K 3 2j{i>j) 



with A = -j= sm9 w , B = ^= cos 6^, and 



m 



V 



■my 



Y U1 .M) Y U11 (H° U ) Y V1 M) Y U13 (H° U )\ 
Y V2i (i>9) Y V21 (H° U ) Y U22 (H° U ) Y V23 (H°) 
Y VZ M) YvM Y, M (H° a ) J 



o \ 


-h(H° u ) 
-XsiH^+YvJvi) 

2K 13j (l>j) 
2K 23j {vfj 

(3.2) 



(3.3) 



For simplicity the summation convention on repeated indices was used in the above 
two equations. The matrix (|3.1| ) is of the seesaw type giving rise to the neutrino 
masses which have to be very small. This is the case since the entries of the matrix 
M are much larger than the ones in the matrix m. Notice in this respect that the 
entries of M are of the order of the electroweak scale while the ones in m are of the 
order of the Dirac masses for the neutrinos |3|, f|] . Therefore in a first approximation 
the effective neutrino mixing mass matrix can be written as 



m 



eff 



—m 



T 



M~ 



m. 



(3.4) 



7 



Because m e // is symmetric and m\^m e ff is Hermitian, one can diagonalize them by 
a unitary transformation 



UuNs m effUMNs = diag{m Vl ,m V2 ,m vl 



(3.5) 



U MNs m lff m effU M Ns = diag(ml 1 ,ml 2 ,ml 3 ). (3.6) 

In Appendix B, Eq. ( B.l| ), we present an approximate analytical expression for 
the effective neutrino mass matrix of the /iz/SSM with SCPV, neglecting all the terms 
containing Y 2 u 2 , Y„v and Y v v z in Eq. (|3.4j ) due to the smallness of Y v and v [Q. In 
the limit of vanishing phases (f Vu = ip Vd = ip u c = ip v . = 0, Eq. ( |B.1|) is reduced to 
Eq. ( |B.9| ). This is the formula that we will use in the following, in order to have a 
qualitative idea of how the seesaw mechanism works in this model. 



Let us first rewrite the expression (B.9) in the following form: 



[ m eff\real)ij 



with 



— 6kX c ^ Y Uj {l 3 Sij , 



2M 



ViVj + 



v d (Y v Vj + Y v .v^ Y V Y V 



3A 



9A 2 



(3.7) 



M, 



eff 



M 



2M (kv c + Xv u v d ) 3Xu c 



2kW 



c 2V u V d . \v* 



+ 



(3.8) 



which coincides with the result in || , where the possibility of obtaining an adequate 
seesaw with diagonal Yukawa couplings was also pointed out. Here v 2 = v \ + v\ + 
^\ v 2 w v 2 + v \ with v ~ 174 GeV has been used, since z/j << v u ,v d 0, and let us 
recall that we are also using couplings Aj = A, a tensor k with terms Km = Ki = k 
and vanishing otherwise, diagonal Yukawa couplings Y m = Y v ., VEVs = u c , and 



m Mi ' ^h^ 



In the limit where gauginos are very heavy and decouple (i.e. M 
Q3.7| ) reduces to 



oo), Eq. 



\JTleff\real)ij 



6 Kf 



^.(1-3^). 



(3.9) 



It is interesting to note that in contrast with the ordinary seesaw (i.e. generated 
only through the mixing between left- and right-handed neutrinos), where the case 
of diagonal Yukawas would give rise to a diagonal mass matrix of the form 



\JTl(.j j ordinary seesaw/ij 



~ V vXvjYvfiij 

2kv c 



(3.10) 



in this case we have an extra contribution given by the first term of Eq. ( |3.9| ). This 
is due to the effective mixing of the right-handed neutrinos and Higgsinos in this 
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limit, and produces off-diagonal entries in the mass matrix. Besides, when right- 
handed neutrinos are also decoupled (i.e. v c — > oo), the neutrino masses are zero as 
corresponds to the case of a seesaw with only Higgsinos. 

Another observation is that, independently on the nature of the lightest neu- 
tralino, Higgsino-like or z/ c -like or even a mixture of them (recall that the v c can be 
interpreted also as the singlino component of the neutralino since R-parity is broken) , 
the form of the effective neutrino mass matrix is the same when the gauginos are 
decoupled, as given by (|3.9|) . 



Another limit which is worth discussing is v c — ► oo. Then, Eq. (|3.7|) reduces to 
the form 



[ m eff\real)ij 



, v d {Y Vi i/j + Y v i/i) Y U Y v 2 ' 
^ + 3A + 



2M 

We can also see that for Vd — > (i.e. tan/3 = ^ — > oo) one obtains 



(3.11) 



{m ef f\real)ij ^ (3-12) 

Note that this result can actually be obtained if V\ » g A , and that this relation 
can be fulfilled with v d ~ v u ~ 174 GeV for suitable values of A. It means that 
decoupling right-handed neutrinos/singlinos and Higgsinos, the seesaw mechanism 
is generated through the mixing of left-handed neutrinos with gauginos. This is 
a characteristic feature of the seesaw in the well-known bilinear R-parity violation 
model (BRpV) @. 

The seesaw in the /zz/SSM comes, in general, from the interplay of the above 
two limits. Namely, the limit where we suppress only certain Higgsino and gaugino 
mixing. Hence, taking Vd — > in Eq. ( |3.7|) , which means quite pure gauginos but 



Higgsinos mixed with right-handed neutrinos, we obtain 



v 2 



(m e ff\real)ij ^ ~ *-Y Vi Y Vj (l ~ 3(%) - ——ViVj , (3.13) 
KsKV Z JWefT 

As above, we remark that actually this result can be obtained ii Vi » g\ . The 
effective mass M e fj = M (l — — ^ e3 J represents the mixing between gauginos and 
Higgsinos- v c that is not completely suppressed in this limit. Expression Q3.13| ) is 
more general than the other two limits studied above. On the other hand, for typical 
values of the parameters involved in the seesaw, M eff w M, and therefore we get a 
simple formula that can be used to understand the seesaw mechanism in this model 
in an qualitative way, that is 

v 2 1 



The simplicity of Eq. (|3.14 ), in contrast with the full formula given by Eq. (|3.7| ), 



comes from the fact that the mixing between gauginos and Higgsinos-z/ c is neglected. 
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AmL/10" 5 eV 2 


sin 2 9 12 


sin 2 9 13 


sin 2 6*23 


AmLJlO" 3 eV 2 


7.14-8.19 


0.263-0.375 


< 0.046 


0.331-0.644 


2.06-2.81 



Table 1: Allowed 3cr ranges for the neutrino masses and mixings as discussed in [27]. 



To continue the discussion of the seesaw in the /ivSSM, let us remind that two 
mass differences and mixing angles have been measured experimentally in the neu- 
trino sector. The allowed 3a ranges for these parameters are shown in Table [I]. We 
also show the compositions of the mass eigenstates in Fig. |TJ for the normal and 
inverted hierarchy cases. For the discussion, hereafter we will use indistinctly the 
subindices (1,2,3) = (e, r). 

Due to the fact that the mass eigenstates have, in a good approximation, the 
same composition of and v T we start considering Y v% = Y V3 and u 2 = ^3, and 




atmospheric 
~2xKT 3 eV 2 



solar~7xlCr 5 eV 2 



solar~7xlO- 5 eV 2 



atmospheric 
~2xl(T 3 eV 2 



Figure 1: The two possible hierarchies of neutrino masses as shown in p8 |. The pattern 
on the left side corresponds to the normal hierarchy and is characterized by one heavy state 
with a very little electron neutrino component, and two almost degenerate light states with a 
mass difference which is the solar mass difference. The pattern on the right side corresponds 
to the inverted hierarchy and is characterized by two almost degenerate heavy states with a 
mass difference that is the solar mass difference, and a light state which has very little electron 
neutrino component. In both cases the mass difference between the heaviest/lightest eigenstate 
and the almost degenerate eigenstates is the atmospheric scale. 
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therefore Eq. ( ft. 14 ) takes the form 




m eff = c A B , (3.15) 



where 



v 2 1 
d = --^Y 2 - —v 2 

3kv c V1 2M u 

vl 1 



§kp cViV2 2M 
v 2 1 

Zkv- ^ 2M 2 ' 
bKv c 2M 



The eigenvalues of this matrix are the following: 



^ (A + B - a/8c 2 + (A + B - d) 2 + d),^(A + B + a/8c 2 + {A + B-d) 2 + d),A-B , 

(3.17) 

and the corresponding eigenvectors (for simplicity are not normalised) are 



A+B+^8c 2 +(A+B-d) 2 -d 

2 ' ' 

-A-B+ v / 8c 2 +(A+B-d) 2 +rf 

(0,-1,1). (3.18) 

We have ordered the eigenvalues in such a way that it is clear how to obtain the 
normal hierarchy for the v Vk -v T degenerate case. Then we see that sin 2 # 13 = and 
sin 2 #23 = \i as in the tri-bimaximal mixing regime. Also we have enough freedom to 
fix the parameters in such a way that the experimental values for the mass differences 
and the remaining angle 9\2 can be reproduced. It is important to mention that the 
above two values of the angles are a consequence of considering the example with v^- 
u T degeneration, and therefore valid even if we use the general formula ( |3.7|) instead 
of the simplified expression ( ft.!4| ). Notice that Eqs. ( ft.l5[) , ( ft,17| ) and ( |3.18| ) would 



be the same but with the corresponding values of A, B, c and d. 

Let us remark that the fact that to obtain the correct neutrino angles is easy 
in this kind of seesaw is due to the following characteristics: R-parity is broken and 
the relevant scale is the electroweak one. In a sense we are giving an answer to the 
question why the mixing angles are so different in the quark and lepton sectors. 
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To show qualitatively how we can obtain an adequate seesaw with diagonal 
neutrino Yukawa couplings, let us first consider the limit 2 c — > . In this limit the 
electron neutrino is the lightest neutrino, and is completely decoupled from the rest. 
The second eigenvector has no v e composition (sin# 12 — > 0), and it is half and 
half v T . Understanding this case we can easy generalized the situation to the case 
sin 9 12 7^ 0, switching on the parameter c. The eigenvalues in this limit are 



d , A + B, A-B, 



(3.19) 



where 



14 = 

\A + B\ 



U u y2 , _J_ 

3kz/ c Vl 2M 



V , 



U y 2 _L J_ 

6kv c V2 M 

2 



1A 



\A-B\ 



v 



(3.20) 



We can see that Am 2 



atm 

2 



\4AB\ 



v i Y i 

M8k 2 i/ c2 

' uj y2 
- 3ki/ c v x 



1 



Y1Mrv c > 



and Am 2 sol 



\(A + B) 2 d 2 \ — ( 6k u uC Y 2 2 , M „ 2J \ 3ki/ c + Ui i 2M' 

It is important to note that we need \A — B\ > \A + B\ for the normal hierarchy 
case, otherwise the 9u angle is zero even when c is not neglected. This is easy to 
obtain for M » 2nv c . If M ~ 2k,u c , using different signs for the effective Majorana 
and gaugino masses helps to fulfill the above inequality. For this to hold with our 
convention, one must take M < 0. 

In the inverted hierarchy scenario \A — B\ > \A + B\ leads the angle 9\2 to 
zero also with c / which is not phenomenologically viable. Then we impose 
\A — B\ < \A + B\. Note that when c is switched on, the parameter d has to be large 
enough for having the associated neutrino with an intermediate mass, as corresponds 
to the inverted hierarchy scenario. Therefore in this case we can also have easily the 
tri-bimaximal mixing regime for M « 2kv c . When M ~ 2kz/ c , having M > helps 
to fulfill the above condition. 

Let us finally remark that we can get the complete tri-bimaximal mixing regime 
sin 2 9i 3 = 0, sin 2 9 23 = 1/2 and sin 2 12 = 1/3 fixing in Eq. (P15|) c = A + B-d. In 
this way we obtain the eigenvalues 



-(A + B) + 2d , 2{A + B) - d, A- B, 
and from Eq. (ft. 18 ), after normalization, we arrive to sin 2 9u = 1/3. 



(3.21) 



2 Actually this limit can be obtained taking Y Vl — * 0, v\ — > 0, implying c — > 0, and also d — * 0, 
and leading to similar conclusions. This limit means that the electron neutrino is decoupled from 
the other two neutrinos, having a negligible mass. 
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Breaking the degeneracy between the Y v and v of the muon and tau neutrinos, it 
is possible to find more general solutions in the normal and inverted hierarchy cases. 
We will show this with numerical examples in the next section, working always in 
the case M ~ 2kv c . Note also that in the case of degenerate v ft -v r parameters, as 
the Dirac CP phase always appears in the MNS matrix in the form sin Oi^e 15 (see eq. 
( f4.2p below), the SCPV effect is suppressed since sin #13 is negligible. This is not the 
case if we break the degeneration between and v T . 

When the vacuum is non CP-conserving the situation is more complicated since 
new relative phases are present, but the idea still holds. In the next section we will 
use the above results to find numerical examples in the general case where also phases 
are generated through complex vacua. Examples where changing the sign of M the 
second and third eigenvalue are interchanged and the behavior is similar to the one 
described in this section. 



4. Results 

In Section we have given a simple argument to show that the /zz/SSM can violate 
CP spontaneously. In Section 3 we have discussed how to obtain correct neutrino 
masses and mixing angles. In this section we sketch the numerical method used for 
the search of global minima of the /iz/SSM with SCPV, giving rise also to an effective 
neutrino mass matrix that reproduces correctly the phenomenology of the neutrino 
sector according to observations. We also give some examples. 

For simplicity, we assume that all the parameters appearing in the potential are 
diagonal in flavor space at the electroweak scale, except the trilinear (A K K)ijk terms 
whose entries different from (A K K)m are relevant to break CP spontaneously. We 
introduce the following notation for the flavor diagonal free parameters of the scalar 
potential: Ki, Y v ,, {A v Y v )i, wif , m h with i = 1,2,3 being flavor indices. Under 
this assumption, the neutral scalar potential in (|2.4| ) is obviously simplified, and as 
a consequence also the minimization conditions ( |A.lf|A.7|) are simplified. In addition 
to the complex VEVs, the potential depends on Aj, Y Vi , (A R K,)ijk, (A\\)i, (A v Y v )i, 
m Hd , m Hu , m„c and m~ L% . 

The strategy followed to find minima of the model consists of solving the mini- 
mization equations in terms of the soft parameters that are linear in those equations. 
More precisely, the three minimization equations ( |A.7|) , corresponding to = 0, 
are used to solve the values of (A v Y v )i. Using this result, Eqs. ( |A.6| ) for i = 2,3, 
corresponding to a 9 Y = 0, are then solved for (A A A) 2 3- Repeating the procedure 

using the equation (|A.5|) , J^- = 0, one obtains (A\X)i. Finally, Eq. (|A.6|) for i — 1 
is used to get (A K K)m. The conditions with respect to the moduli of the VEVs 
( |A.1[ - [A.4D are used to get the squared soft masses. Once this is done, we ensure that 



the critical point found (i.e. with non-vanishing phases for the VEVs) is a global 
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minimum through a numerical procedure. As discussed in one has to check in 
particular that the minimum found is deeper than the local minima with some or all 
the VEVs vanishing. 

To accomplish the numerical task of finding global minima we need as inputs 
the eight moduli and seven phases of the VEVs, the Aj, Ki and Y v . couplings and 
the soft-trilinear terms (A K K)ij^ with (i,j,k) ^ (1,1,1). For simplicity, we assume 
a special structure for the latter: (v4 k k) 2 22 = {,A k k)^z, a common value for (A^k)^ 
with i,j, k 7^ 1, and another common value for {A K K)ijf. with one or two indices 
equal to 1. Moreover, let us recall that the modulus of the SUSY Higgs VEVs, can 
be determined from v 2 = v\ + v 2 + ^ uf « v\ + v\ with v ~ 174 GeV, and the value 
of tan/3 = 2*. 

Once we find global minima, the next step is to build the neutralino mass matrix 
and to diagonalize it perturbatively in order to extract the effective neutrino mass 
matrix. Diagonalizing the effective neutrino mass matrix, we can extract the mass 
differences and the mixing angles of the neutrino sector and compare them with the 
data. The key for obtaining a phenomenologically viable neutrino sector, once we 
are in a global minimum, consists of varying either the neutrino Yukawa couplings, 
the left-handed sneutrino VEVs or the soft gaugino masses. This approach does not 
alter the vacuum structure previously obtained. 

Let us now describe the details on how we proceed with the phenomenological 
analysis of the neutrino sector of the model. First, we assume for simplicity the GUT 

2 

inspired relation between the gaugino masses M\ and M 2 , M\ = -7M2, implying 
Mi ~ 2 Mi at low energy. As discussed in Section 3, one has to diagonalize the 
neutrino effective mass matrix, m e ff = —m T -M~ 1 -m. Since it is a complex symmetric 
matrix, it can be diagonalized with an unitary transformation, as it is shown in Eqs. 



75]) and ( p.6|) . For the MNS matrix we follow the standard parameterization 

U MNS = diag(e i5 % e w ", ■ V ■ diag(e" i ^/ 2 , e"^ 2 , 1), (4.1) 
where 4>\ and 2 are the Majorana phases and V is given by 

Ci 2 ci3 S12C13 si 3 e~'' 

V = ( — C23S12 — S23Si3Ci2e ,<5 C23C12 — S23Si3Si2e ,<5 S23C13 | . (4-2) 
S23 s 12 — C23Si3Ci2e* 5 — S23C12 — C23Si3Si2e i<5 C23C13 

Here = cos 8^ and = sin#y whereas 5 is the Dirac CP violating phase. The 
conventions used for extracting the mixing angles and the Majorana and Dirac phases 



from Eqs. ( f4.1|) and ( |4.2j ) are outlined in Ref. [p9|j . 

Taking all the above into account, we show in Table |2| the parameters that 
characterize an example of a global minimum that breaks CP spontaneously. The 
values of the soft parameters not determined by the minimization equations have 
been chosen to be (A K K)m = 280 GeV for i ^ 1, (A K K)ijk = —40 GeV for i, j, k ^ 1, 
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Xi = 0.13 


= 0.55 


i/? = 1000 GeV 


tan (3 = 29 


ip v = -TT 


<P*t = f 


^ = tp** = -f 


xi = -f 


X2 = X3 = f 



Table 2: Numerical values of the relevant input parameters for a global minimum that 
breaks CP spontaneously. 



Y U1 = 4.25 x 10~ 7 


Y U2 = Y vz = 1.36 x 10~ 6 


M 1 = -340 GeV 


v x = 3.88 x 10~ 5 GeV 


v 2 = u 3 = 1.24 x 10~ 4 GeV 





Table 3: Numerical values of the neutrino/neutralino inputs that reproduce the neutrino 
experimental constraints, and correspond to the normal hierarchy scenario. 



{A V Y V ) X ~ -0.0031 GeV 


(A v y v ) 2 ~ -0.010 GeV 


(A U Y U ) 3 ~ -0.010 GeV 


(A X X) 1 ~ -1487 GeV 


(A A A) 2 ~ -679 GeV 


(A A A) 3 ~ -679 GeV 


{A K K) in ~ -0.25 GeV 


m 2 Hd ~ 7.0325 x 10 7 GeV 2 


m 2 ff ~ -47200 GeV 2 


m\ ~ 260140 GeV 2 


m 2 ~ -100820 GeV 2 


mlc ~ -100820 GeV 2 

"a 


m 2 ~ m 2 = m 2 = 137 x 1Q 7 Ge y2 

Li | l^a 





Table 4: Values of the soft terms calculated with the minimization equations for the global 
minimum associated to the parameters shown in Table |2L 



and (A K K,)ijk = —120 GeV for one or two indices equal to 1. In Table |3| we show 
the neutrino/neutralino inputs used in order to obtain a i/ M -v T degenerated case with 
normal hierarchy, producing values of masses and angles within the ranges of Table 
[I]. In particular, we obtain sin 2 6*13 ~ and sin 2 #23 = 0.5, as expected from the 
discussion in Section 3, sin 2 6 12 = 0.323, and neutrino masses mi = 0.00305 eV, 
m 2 = 0.00949 eV and m 3 = 0.05091 eV, producing Am 2 solar = 8.08 x 10~ 5 eV 2 and 
Am 2 tm = 2.50 x 10~ 3 eV 2 . The corresponding values of the soft terms calculated 
with the minimization equations are presented in Table |j. 

It is worth noticing that for this solution, the soft masses of the left-handed 
sneutrinos, m^., do not need to be very different, and, actually, in this case they are 
almost degenerate ~ 3700 GeV. This can be understood using the minimization equa- 

Y Y u 

tions ( |A.4| ), neglecting the terms with products of Yukawas. When — = — , V i,j, 
one obtains m? = m| . However, we have to point out that the values obtained for 
other soft parameters are not so natural in a SUSY framework. Notice for example 
that A u ~ —7 TeV, A Xl ~ — 11 TeV, whereas A K111 ~ —0.5 GeV. Indeed, this is a 
consequence of the particular solution shown in Table @. 

Although it is non-trivial to find realistic solutions, since many minima which 
apparently are acceptable, at the end of the day turn out to be false minima, we 
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Xi = 0.10 


«i = 0.35 


!/j = 835 GeV , ^ 2 C = i/g = 685 Gel/ 


tan p = 29 


= - 71 " 


</M = f 


<p*s = ^ = -f 


Xi = "f 


X2 = X3 = f 



Table 5: Numerical values of the relevant inputs for the second global minimum discussed 
in the text, that breaks CP spontaneously. 



Y Ul = 5.4 x 10" 7 


Y U2 = Y vz = 9.2 x 10~ 7 


M x = -340 GeV 


i/i = 3.7 x 10- 5 GeV 


v 2 = v 3 = 8.8 x 10~ 5 GeV 





Table 6: Numerical values of the neutrino/neutralino inputs for the second global min- 
imum discussed in the text, that reproduce the neutrino experimental constraints and 
correspond to the normal hierarchy scenario. 



(A u Y u )i ~ -0.00209 GeV 


(A V Y V ) 2 ~ -0.00294 GeV 


(A U Y U ) 3 ~ -0.00294 GeV 


(A\\)i ~ -156 GeV 


(4 A A) 2 ~ -84 GeV 


(A A A) 3 - -84 GeV 


{A k k) x11 ~ 12.7 GeV 


mf^ ~ 5.36 x 10 6 GeV 2 


mi ~ -37910 GeV 2 


m| f ~ 51035 GeV 2 


m? § ~ 69155 GeV 2 


m% ~ 69155 GeV 2 


m 2 = 8.07 x 10 6 GeV 2 


m 2 ? = 3.92 x 10 6 GeV 2 


m 2 ? = 3.92 x 10 6 GeV 2 



Table 7: Values of the soft terms calculated with the minimization equations for the second 
global minimum discussed in the text, associated to the parameters shown in Table |5[ 



have been able to find more sensible solutions. This is the case of the one shown in 
Table [5|, with the values of the input soft parameters {A^m = — 150 GeV for % ^ 1, 
(A K K)ijf. = 75 GeV for k ^ 1 and (A K K)ij k = —50 GeV for one or two indices 
equal to 1. For example, lowering the values of v c one is able to lower the trilinear 
terms A v ~ —3 TeV in order to fulfill Eqs. ( |A..7Q (also lowering k contributes to this 
result), and also the soft masses m^. ~ 2.8 TeV, as shown in Table [7]. Lowering A 
one is able to lower the trilinears A\ 1 ~ —1.5 TeV, Ax 23 ~ —840 GeV, in order to 
fulfill Eqs. (|A.5| ) and (|A.6| ). Notice finally that the use of non-degenerate v\ allows 
to increase the trilinear A K111 ~ 36 GeV. In Table |^ we show the corresponding 
neutrino/neutralino inputs producing values of masses and angles within the ranges 
of Table 0. 

Modifying the values of the angles we can also obtain other interesting solutions. 
See for example the one shown in Tables [8], || and 0. In this case the values of the 
input soft parameters are chosen to be {A K K)in = —200 GeV for i ^ 1, (A^k)^ = 125 
GeV for i,j,k ^ 1 and (A K K,)ijk = —75 GeV for one or two indices equal to 1. 
Notice that now the values obtained for the soft terms are also of this order. In 
particular, the trilinears are A Vl ~ —657 GeV, A U2 „ ~ —429 GeV, A\ 1 ~ —990 GeV, 
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Xi = 0.10 


Ki = 0.42 


i/« = 850 GeV , i$ = i% = 550 Gel/ 


tan (3 = 29 


V9 t , = -7T 


¥M = f 


<Pu* = (P^ = "f 




X2 = X3 = f 



Table 8: Numerical values of the relevant inputs for the third global minimum discussed 
in the text, that breaks CP spontaneously. 



Y Ul = 1.9 x 10" 7 


Y U2 = Y V3 = 8.5 x 10~ 7 


M x = -100 GeV 


v x = 6 x 10~ 5 GeV 


U2 = v 3 = 4.9 x 10~ 5 GeV 





Table 9: Numerical values of the neutrino/neutralino inputs for the third global minimum 
discussed in the text, that reproduce the neutrino experimental constraints and correspond 
to the normal hierarchy scenario. 



(A V Y V ) X ~ -0.000125 GeV 


{A U Y U ) 2 ~ -0.000365 GeV 


(A U Y U ) 3 ~ -0.000365 GeV 


{A\X)i ~ -99 GeV 


(^ A A) 2 ~ -83 GeV 


(^ A A) 3 ~ -83 GeV 


{A k k) 11x ~ 41.9 GeV 


m\ d ~ 3.6 x 10 6 GeV 2 


m 2 . ~ -25118 GeV 2 

n u 


m%o ~ -24393 GeV 2 


m| § ~ 208377 GeV 2 


m 2 -c ~ 208377 GeV 2 


m 2 f = 394777 GeV 2 


m? = 903528 GeV 2 

Li 


m 2 f = 903528 GeV 2 



Table 10: Values of the soft terms calculated with the minimization equations for the third 
global minimum discussed in the text, associated to the parameters shown in Table |8[ 



A\ 2 3 ~ —830 GeV, and A Klll ~ 100 GeV. For the soft masses we obtain ~ 628 
GeV, m~ L2 3 ~ 950 GeV. 

A general analysis of the parameter space, finding other interesting complex 
vacua, is obviously extremely complicated given the large number of parameters in- 
volved, and beyond the scope of this paper. Nevertheless, we have checked that other 
sensible solutions can indeed be obtained modifying adequately the parameters. In 
the following we will work with the solution associated to the parameters of Table 0, 
since the discussion below is essentially valid for other solutions. Our strategy will 
consist of varying the neutrino/neutralino inputs Y v ., and M\ in such a way that 
the derived neutrino mass differences and mixing angles are within the ranges of 
Table 1. As mentioned above, this procedure will not alter the vacuum structure 
found. Notice in this respect that gaugino masses do not contribute to the mini- 
mization equations, and that the values of Y u . and z/j are very small. Let us also 
mention that this strategy can indeed be applied to the much more simple issue of 
analyzing real vacua. In particular, it was shown in Q that many global minima 
with real VEVs can be found. For them neutrino/neutralino inputs Y Vi , Ui, M x , 
similar to those studied here are also valid. 
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As noted in Sect. |3| we have chosen M\ < in order to guaranty a viable 9 12 
angle. It is worth pointing out here that a redefinition of the parameters leaving the 
Lagrangian invariant can be made, in such a way that Mi becomes positive and other 
parameters such as the VEVs become negative, describing indeed the same physics. 
In our convention the VEVs, vj, v u , z/ c , are always taken positive. 

We would also like to stress that all the numerical results have been obtained 
without any approximation, that is, with the exact expression of the 10 x 10 neu- 
tralino mass matrix, calculating numerically the effective neutrino mass matrix and 
diagonalizing it. 

Let us first study how the neutrino mass differences depend on the inputs. In 
Sect. |3] we showed that in this scenario there are two different contributions to the 
seesaw mechanism; the one involving right-handed neutrinos (and Higgsinos) given 
by ^ " ; where the Dirac and Majorana masses are parameterized by Y u .v u and 
2ku°, respectively, and the contribution coming from the gaugino seesaw given by 
+ ^2 ' wnere the Dirac and Majorana masses are parameterized by g a Vi and 
M a , respectively, with a = 1,2. 

Figs. Qa and 0b show that the heaviest eigenvalue (dashed line) has very little 
electron-neutrino component, as expected in the normal hierarchy scenario (see Fig. 
[I]), and therefore it does not depend on (Y Vl v u ) 2 /(2kv c ), whereas the intermediate 
(solid line) and lightest (dotted line) eigenvalues, that have sizeable electron-neutrino 
components, grow with this term. As a consequence of the latter, the squared solar 
mass difference grows as well. On the other hand, following the arguments related to 
Eq. (|3.14j ), we can see in Figs. ||c and |^d that the heaviest eigenvalue is controlled 
by the contribution of the seesaw with right-handed neutrinos having an important 
muon/tau neutrino composition, thus we observe how the heaviest eigenvalue grows 
with (Y U2 v u ) 2 / (2kv c ) and, as a consequence, the squared atmospheric mass difference 
grows accordingly. The variation with (Y V3 v u ) 2 / \2nv c ) is analogue. 

Fig. I is analogous to Fig. || but showing the squared neutrino mass differences 
dependence on the gaugino seesaw component. In this case, because the heaviest 
eigenstate (dashed line) practically does not mix with the electron neutrino we can 
see that it does not vary with ((giis^/Mi + (g 2 Pi) 2 / M 2 ) 2 for i = 1,2,3. On the 
other hand, the intermediate eigenstate grows with the mixing with the gauginos, as 
explained in Sect. |3| with Mi < 0, therefore the squared solar mass difference also 
grows. 

Let us now discuss the mixing angles. Note that in the v il -v T degenerate case with 
normal hierarchy and M\ < we have obtained sin 2 #13 = and sin 2 #23 = §• In Fig. 
^ we present the variation of sin 2 9yi with the ratio of the parameters that control 
the gaugino seesaw, fr 2 /& 2 , where for the sake of simplicity we take hi = Y Ui Vd + 3A^ 
and we do not consider the complicated factors containing phases in Eqs. ( |B.5| ). 

To obtain results different from sin 2 6*23 ~ | and sin 2 #13 ~ 0, in the following 
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Figure 2: Squared neutrino masses versus (Y v .v u ) / (2kv c ) 2 . (a) and (b) show for i = 1 the 
two heaviest and lightest neutrinos, respectively. The same for (c) and (d) but for i = 2. 



we consider the possibility of having different values for the Y v and v parameters 
for fi and r neutrinos. We show in Fig. [5]a sin 2 623 as a function of the ratio of 
the term that controls the Higgsino-^ c seesaw, a^/a 2 . When a M /a T goes to 1, the 
v^-Vt degeneracy is recovered and sin 2 #23 goes to 1/2 as expected. In Fig. |b we 
show sin 2 613 as a function of 4aM " T „ that is a good measure of the degeneration in 
this case. Note that when 4a /i a T /(a^ + a T ) 2 



sin 2 fl 
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1 the degeneracy is recovered and 
as expected. The parameters aj have been defined in Eq. ( B.2 ). Let us 
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Figure 3: The same as in Fig. 2 but for the squared neutrino masses versus [(giVi) 2 /Mi + 
(g2Vi) 2 /M 2 }\ 



point out that sin 2 6>i 3 < 1CT 3 since we are breaking the degeneration between /i and 
t neutrinos but the term that controls the higgsino-z/ c seesaw for the first family is 
very small compared to the other two families. 

As mentioned previously, the /x^SSM with SCPV also predicts non-zero CP 
phases in the MNS matrix. We have checked numerically that for each of the ex- 
perimentally allowed regions found, the two Majorana CP phases and the Dirac CP 
phase are different from zero. This fact is reflected in Fig. ^| where we present two 
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Figure 5: (a) The variation of sin 2 6*23 with respect to the relevant term that controls its 
evolution, a 2 /a 2 . (b) The variation of sin 2 #13 with respect to the term that measures the 
v v -v T degeneracy. 



plots in the 5 — <fi\ and 5 — 02 planes (Dirac-Majorana CP phases) constructed vary- 
ing all the inputs in the neutrino sector. However, it is fair to say that due to the 
smallness of sin 2 13 ~ 1CT 3 in this region, the CP violation effects of the phases of 
the VEVs turn out to be suppressed in the MNS matrix because the Dirac CP phase 
always appears in the form sm9i 3 e tS . 

In order to complete the discussion about the neutrino sector in this scenario, 
we will consider the possibility M 1 > instead of Mi < . In Sect. 3 we have seen 
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Figure 6: 6 — <pi plane (a) and 5 — 4>2 plane (b) for the scenario with normal hierarchy and 
negative gaugino masses M < 0, varying simultaneously Y v .,Vi,M\. 



that with Mi > it is more complicated to have a degeneracy between muon and 
tau neutrinos because it is easy to obtain sin 2 #12 ~ 0, in contradiction with the data 
(see Table |Tj). Thus we will show a region where breaking the degeneracy z/ M -z/ T a 
normal hierarchy is obtained with M\ > 0. This region is around the point of the 
parameter space shown in Table [XT . In this example the angle sin 2 # 13 can easily be 
made small as required by the data, but it is not necessarily negligible. Thus the CP 
violating effects would be present in the MNS matrix. Besides, we can roughly say 
that sin 2 #13 and sin 2 6*12 are interchanged with respect to the case discussed above 
with Mi < 0. For completeness, in Fig. 0a we show the variation of sin 2 6*13 with 
respect to the term that controls the gaugino seesaw relevant in this case, namely 
6 2 / (6 2 + 6 2 ). We also plot in Fig. [7|b sin 2 Q 12 as a function of the relevant term that 
controls the Higgsino-i^ seesaw ° T As mention above, an interesting feature 
of this region of the parameter space is that the effect of the Dirac CP phase in the 
MNS is not removed, since the value of sin #13 is not negligible. Fig. [8] shows the 
derived CP phases of the MNS matrix. 



Y U1 = 9.54 x 1(T 7 


Y V2 = 9.47 x 10- 7 


Y U3 = 2.31 x 1(T 7 


Mi = 350 GeV 


v x = 8.59 x KT 5 GeV 


u 2 = 2.25 x 1(T 4 GeV 


z/ 3 = 2.29 x KT 4 GeV 





Table 11: Numerical values of the relevant neutrino/neutralino-sector inputs that re- 
produce the neutrino experimental constraints, and correspond to the normal hierarchy 
scenario with Mi > 0. 
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4 a T /(a^+a T ) 

(a) (b) 

(a) The variation of sin 2 #13 with respect to the relevant term that controls its 
(b) The variation of sin 2 #12 with respect to the relevant term 4a fJ ,a T /(a IJ , + a T ) 2 . 



For the sake of completeness, we show in Table [12] an example where the inverse 
hierarchy scenario is achieved. 

At this point it is clear that there are many regions with different characteristics, 
different compositions for the lightest neutralino or regions close to the tri-bimaximal 
mixing regime for normal or inverted hierarchy that can be found with different neu- 
trino parameters. Furthermore, we have seen that the /i^SSM with SCPV predicts 
non-zero CP-violating phases in the neutrino sector. If in the future a non-zero CP 
violating phase in the lepton sector is measured, SCPV as the one analyzed here 
could be a possible source. 

Neutrino oscillations are sensitive only to the Dirac CP phase (insensitive to the 
Majorana phases). Let us briefly comment about the possible determination of 5 in 
future neutrino experiments. The conservation of CP implies P[y a — > v$) = P(v a — > 
Dp). If CP is not conserved, we would have [|30[ 



P(v» -> Ve) ~ P{pn -> *«) = -16Jsin ( — ^- ] sin ( — ^_ ] sin ' 



AE \ AE \ AE 



(4.3) 



where L is the oscillation length, E is the neutrino beam energy and J is the Jarlskog 
invariant for the neutrino mass matrix which is given by J = Si2Ci2S23C23Si3C 2 3 sinrj. 
There is only an upper experimental limit for J, J < 0.04. The reason is that J 
depends on #13 and 5, which are currently unknown. If 813 vanishes (recall the bound 
sin 2 #13 < 0.038) J vanishes and the effect of CP violation via ( f4.3|) would be unob- 
servable. The same occurs if there was a degeneracy in the neutrino masses. In spite 
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Figure 8: 5 — <p\ plane (a) and 5 — 4>2 plane (b) for the scenario with normal hierarchy and 
positive gaugino masses M > 0, varying simultaneously Y v .,Ui, M\. 



Y Vl = 5.98 x 1CT 7 


Y U2 = 1.32 x irr 6 


Y U3 = 1.40 x 10" 6 


M x = 340 GeV 


ui = 3.276 x 10~ 4 GeV 


v 2 = 6.20 x lCT 5 GeV 


u 3 = 6.56 x 10" 5 GeV 





Table 12: Numerical values of the relevant neutrino/neutralino inputs that reproduce the 
neutrino experimental constraints, and correspond to the inverted hierarchy scenario. 



of these extreme situations the process ( |4.3j ) implies that long baseline experiments 
allow the observation of CP violation due to the Dirac phase S in the neutrino sector. 
Two experiments are designed for this purpose: NOz/A plj and the T2KK detector 



On the other hand, although Majorana phases affect neutrinoless double beta 
decay 0vf3f3 |53|], their determination turn out to be difficult. 



Let us finally briefly discuss the implications of the CP-violating phases con- 
cerning the electric dipole moments (EDMs) in the /^z/SSM. As is well known, EDMs 
impose important constraints on supersymmetric theories. The MSSM (with explicit 
CP violation in the soft Lagrangian) predicts EDMs about three orders of magnitude 
larger than the experimental bounds for the EDM of the electron and neutron if the 
supersymmetric CP violating phases are 0(1) and the supersymmetric particles have 
masses near their current experimental bounds (9(100 GeV) ]34]. There are three 



kind of solutions to this problem in supersymmetric theories. First, if the super- 
symmetric CP violating phases are very small, of order O(10~ 2 — 10 -3 ) the EDM 
bounds can be easily satisfied [Bl] . Second, if the supersymmetric scalar particles are 
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Xi = 0.13 


«j = 0.55 


v\ = 900 GeF , z/ 2 c = i/§ = 600 GeV 


tan (3 = 29 


= o 


= m 


= ^ = -IM 


^ ~~ _ 90 


7Y 

X2 — X3 — go 



Table 13: Numerical values of the relevant inputs of a global minimum that breaks CP 
spontaneously with small phases. 



Y Vl = 1.9 x 10" 7 


Y U2 = Y vz = 1.06 x 10" 6 


M x = 300 GeV 


vx = 1.54 x 10~ 4 GeV 


u 2 = u 3 = 2.4 x 10" 5 GeV 





Table 14: Numerical values of the neutrino/neutralino inputs that reproduce the neutrino 
experimental constraints for the global minimum with small phases. 



decoupled with masses larger than about 3 TeV, and thus out of reach of the LHC, 
but not spoiling the solution of supersymmetry to the hierarchy problem, the EDM 



bounds could also be accomplished p5| . Third, there can be internal cancellations 



between the different contributions to the EDMs I0B1 . 



We would like to point out that the /zz/SSM with SCPV could implement these 
three kind of solutions. First of all, the possibility of small supersymmetric CP 
phases is present in our model. Let us show for example a global minimum that 
break CP spontaneously with O(10~ 2 ) CP phases (we have also found global minima 
with 0(1O~ 3 ) phases). The values of the soft parameters not determined by the 
minimization equations are chosen to be (A k k)m = —175 GeV for i ^ 1, (A K K)ijk = 
100 GeV for k ^ 1 and (A K n)ij k = —100 GeV for one or two indices equal to 1. 
The numerical values of the phases and the rest of input parameters are presented 
in Tables O and IL1 . 



It is worth remarking here that in models with SCPV small phases are not 
unnatural, since they arise as a consequence of the minimization conditions (notice 
that the use of phases as inputs in this work is just an artifact of the computation) 
for particular values of the soft terms. 

The other two solutions, heavy scalars and internal cancellations, can also be 
implemented. Notice that the following soft parameters remain free in our model 
because they are not included in either the neutral scalar potential or the neutrino 
sector: (A u Y u )ij , mL , M 3 , (A e Y e )ij , m% . Thus, the solution with heavy scalars 

ij ij 

remains valid for scalar masses heavier than about 3 TeV. We also expect the internal 
cancellation solution to be valid in our model. This is because these free parameters 
enter in the calculation of the EDMs, and we will have enough freedom to find values 
where such cancellations can be accomplished, fulfilling the EDMs bounds. 
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5. Conclusions 



In this work we have studied in detail the neutrino sector of the /xi/SSM. We have also 
shown that, even if all parameters in the scalar potential are real, SCPV is possible 
at tree level, and we have used these vacua to show how a complex MNS matrix can 
arise. 

In particular, we have calculated first the scalar potential of the /zz/SSM with 
real parameters, assuming the most general situation where the VEVs of Higgses and 
sneutrinos can be complex. We have shown, using a simple argument, that CP can 
actually be spontaneously violated in this model. 

Then we have discussed the neutralino-neutrino mass matrix of the /xz/SSM, and 
we have shown how to obtain from it the effective neutrino mass matrix. Although the 
discussion is general, we have applied it also to the particularly interesting case of real 
vacua. We have analyzed how the electroweak seesaw mechanism works in the /iz/SSM 
using approximate analytical equations, particularized for certain interesting limits 
that clarify the neutrino-sector behavior of the model. In addition, we have given 
the qualitative idea of how to find regions in the parameter space of the model that 
satisfy the neutrino experimental constraints. Let us remark that these constraints 
can be fulfilled even with a diagonal neutrino Yukawa matrix, since this seesaw 
does not involve only the right-handed neutrinos but also the MSSM neutralinos. 
Actually, to obtain the correct neutrino angles turns out to be easy due to the 
following characteristics of this seesaw: R-parity is broken and the relevant scale is 
the electroweak one. In a sense, this gives an answer to the question why the mixing 
angles are so different in the quark and lepton sectors. 

Finally, we have presented our results describing the method to obtain numeri- 
cally global minima with SCPV, and giving examples of such minima. Let us empha- 
sized however that, unlike the case with real VEVs where many global minima can 
be found, for the case with complex VEVs such minima are not so easy to find. In 
particular, one has to choose carefully the parameters of the model. For the exam- 
ples found we have shown the dependence of the neutrino mass differences (for both 
normal and inverted hierarchies), mixing angles, and CP phases of the MNS matrix, 
in terms of the relevant neutrino inputs. Last but not least, we have checked that 
different regions of the parameter space can reproduce the neutrino experimental 
constraints. In this context, future neutrino experiments could be able to measure a 
non-zero Dirac CP-violating phase, opening the possibility to SCPV in the /zz/SSM 
as the dominant source. 
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A. Minimization Equations 



Here we write first the eight minimization conditions with respect to the moduli Vd, 



I Y ViVi +V d- V l) V d + m H d V d + V d vl J^(Ai) 2 - ^2(A X \)iU?V u C0S(V?„ + Ifiuf) 
\ i / i i 

+ ^2 V dW V< i V ] COs(<Pv? - Vv-) ~ ^ K ikj^kV U ViVj C0S(^C + (p„c - tp v ) 

~ Y Y ^M v ^ v ] v l cos(xi + (p^ - (fv* - (p v ) ~ Y Y Y "n X i v l u i cos (^ - Xi) = °> 

i,j,k i j 

(A.l) 



- ^ [Y ViVi + V d ~ Vu + m Hu V u + V uV 2 d Y^ 

+ ^2(A V Y V ) 

+ ^ K^jVuV C i v< j COs(<Pvf ~ <Pvj) - Y Y Kijk^kVdtftf COS^c + ip u c - if v ) 

+ YY Y ^i K nkVjV C i v l cos(^c + ip u c - X j) ~YY ZYvij^jVdVuVi cos((p v - Xi) 

i,j,k I i j 

+ Y Y Y »ik Y Vjk v uViVj cos(xi - Xj) + Y Y Y ^y^ v uv1v c j cos(y^ f - (Pvj) = 0, 

i,j k i,j k 

(A.2) 
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^2mlcUjCos(ip u c - (p u c) - (A x X)iV u v d cos((p v + ip u c) + ^(A l/ Y l/ ) ji v j v u cos(xj + ¥„?) 

3 3 

3,k j 

+ Y X ^ X 3 U j V l C0S(^C - Vu c) + ^ Y 2K irnkKlmjtykVl COs(^V? + ip v c - (p„c - ip„c) 
3 j,k,l m 

~ Y Y 2K ijkhVdV u ^ COs(if u c + cp v c — <^«) + 2Y "jl K iM V uVjVk COs(^ + (f u c - Xj) 

3 k j,k I 

~ Y Y y 3 i X kVdVjV C k COs(x i + (fi„c - (p„c -(p v )-^2 Y Vkj \iV d v k v c - cos(xfc + Vv- - <Pu9 - <Pv) 

j,k j,k 

+ Y Y ^y^ v 3 v i v i c ° s (xj -Xk + <pv? - <Pvf) + Y Y Y »ki Y »kA u j cos (^< - <pv?) = o, 

j,k,l j k 

(A.3) 



3 3 3 

+ Y Y Y va K 3lkVuV]v C k C0s(</?„c + ^ - Xi) 
j",fc I 

~ Y Y ^3 X kVdV]Vk COs(Xi + Vv] - ^vl -lf v )-Y Y ^ X 3 V dvl COs(lf v - Xi) 
j,k 3 

+ Y Y i3 Y »kl"jVktf COs(Xi - Xk + ^ ~ (fuf) + Y Y Y ^ Y v^ V l V 3 COS (Xi - Xj) 
j,k,l j k 

= 0. (A.4) 



The seven minimization conditions with respect to the phases <p v , ip u c and x% are: 



- YY ^ijkhvdVu^i^ sin(^c + (p u c - tp v ) 

i,j k 

~ 2 E Y »ii X k v dViVjVk sin (Xi + ^ - ¥>*g - y "y V^z/; sin(v?„ - Xi)] 
+ 2 E(^AA)i^fudU« sm(^„ + ip v c) = 0, (A.5) 
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j 

~ Y X i X i v d^ u j sin (Vuf ~ (Pv°) ~ Y X i X J v l u i u j sinO„! - ^9) 
j j 

- 2 ^ ^ Y^KukVuVjvlvl sin(^c + - Xi ) 

+ ^ ^A^j^X sin(xj + ¥V? - - y«) - ^ Y Vk \iV d Vkv]vt sinfofc + V 2 ^ - - 
~ Y Y ^ Y ^i u i v i u k^i sin(xj - Xfc + </V - Wf) - ^ Y Y »ui Y v kj v l v i vC j sin (</V - ¥v) 

j,k,l j k 

+ (A x X)i^v d v u sm(Lp v + ip y c) - ^(v4 K K) ijA .z/ J c z/jz/^sin(v?^ + (p^ + (p v c) 

- Y( A ^)ji v ^ u j u i sin (Xj + ¥u<r) = 0, (A. 6) 

j 

j 

+ Y Y Y m K 3ikV u ViVjvZ sm(<f u c + y v% - Xi ) + Y Y m X k^i^i sin(x i + <Pvs - <Pv° - <p v ) 

j,k I j,k 

~ Y Y ^ 3 X i v dV 2 u Vi sin(p w - Xi) 
j 

- Y Y ^ Y ^ v ] v ^ v i sin (^ - Xk + ^ - (fpf) 

- Y Y Y ^ Y ^u v l v i u 3 sin (^ - Xj) 

j k 

- y^ j {A v Y v ) ij v u v i v C j sin(xi + y>„°) = 0. (A.7) 



B. Analitical formula for neutrino masses 



The formula presented here is obtained from Eq. (|3.4j ) neglecting terms proportional 
to Y 2 v 2 , Y^v and Y u u 3 , and has been particularized for the simplified case discussed 
in Sect. 4 of a common value of couplings Aj = A, a tensor k with terms Km = = k 
and vanishing otherwise, diagonal Yukawa couplings Y u .. = Y u ., and a common value 
of the VEVs uf = v c . The phase structure of the global minimum discussed in 
Section 4 for analyzing the neutrino sector, <p v c = —ip u c = —(p v c = <p u c and <p L 



29 



— ip V2 = —ip V3 = <fu, has also been used in the computation. Then we arrive to the 
following formula: 



where the parameters have been defined as 

A = {e ilp » c + 2e i3 ^ c )A> 2 + v 2 d f + (8e^ c + & 3lp » c )\Kv c2 v d v u e~ 1 ^ 

- (16 + 16e i2 ^ c + 4e i4 ^" c )MA 2 /tz/ c3 - (8 + 20e i2t ^ c + 8e i4lp " c )M\ 3 v c v d v u e ilpv , 
Z = e^ c [-4e^ c (2 + e i2 ^ c )/tz/ c2 ^ u + e^*A(4M(2 + e i2t ^ c ) 2 /tz/ c3 

- e^- c (1 + 2e l2 ^ c )(^ + w 2 ) 2 ) + 4e l2( ^ c+ ^ ) A 2 Mz/ c ^ u (5 + 4 cos 2</v0], (B.2) 

w "h i = ^ + 4' 

T n = 2e i2 ^[-4e i2(v " c+ ^ ) (2 + e ,2 ^ c )MAV c Wfi + 4e^- c /tz/ c2 ^ u 

+ e^A(-4(2 + e l6 ^ c )MKu c3 + e i3 ^ c (v 2 u + v 2 d )% 
T 22 = T 33 = 2e i(v " c+2v " ) [_4e i2( ^ c+ ^ } (2 + e i2 ^ c )MAV c uA + 4e i3 ^ c /tz/ c2 ^ u 

+ e^A(-4(l + e l2 ^ c + e i4 ^ c )M/tz/ c3 + e i3 ^ c {v 2 u + v 2 ) 2 )], 
T 12 = T 13 = - e i2 ^[-4e i2( ^ c+ ^ ) (2 + e i2 ^ c )MX 2 u c v d v u + Ae i3 ^ c ku° 2 v d v u cos(2yv) 

+ e^ c+ ^)A(4(-3cos(3^ c ) +^in(3^))^^ c3 + + ^) 2 )L 
T 23 = -e i2 ^[-4e i2(2 ^ c+v " ) (2 + e i2 ^ c )MX 2 u c v d v u + 4e i3 ^ c Kv c \ d v u 

+ e^A(-4(-l + 4e l3 ^ c cos(yv))M/tz/ c3 + e l5 ^> 2 + v 2 ) 2 )], (B.3) 

and 

X u = 2Ku c \b 11 ) 2 + 2\u c v d v u e^(b' n ) 2 + e n , 

X 22 = 2Kv c \b 22 f + 2Xu c v d v u e l ^(b' 22 ) 2 + e 22 , 

X 33 = 2Ku c \b„) 2 + 2Az/ c ^ u e^(6 33 ) 2 + e 33 , 

X 12 = 2k I / 3 (& 11 )(& 22 ) + 2\v c v d v u e^(b' 12 ) 2 + e 12 , 

X 13 = 2Ku c \b n )(b 33 ) + 2Az/ c ^ u e^(6' 13 ) 2 + e 13 , 

X 23 = 2/tz/ c3 (6 22 )(6 33 ) + 2Az/ c ^ u e^(6 23 ) 2 + e 23 , (B.4) 
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with 



(6n) = (2 + e i2 ^)\e-^v 1 + e i2 ^v d Y Vl , 
(622) = (2 + e l2 ^)\e l ^u 2 + v d Y U2 , 
(633) = (2 + e l2 ^)\e^v 3 + v d Y 1/ . i , 
(b' u ) 2 = (2 + 5e i2 ^ c + 2e i4 ^ c )AV i2 ^z/ 2 

+ (2 + 2e i2 ^ c + 2e A ^ c )\v d e- iLpv v x Y Vl + e^u^, 
(b' 22 ) 2 = (2 + 5e l2 ^ c + 2e l4 ^ c )AV 2 ^z/ 2 

+ (1 + 4e i2t ^ c + e i ^ c )\v d e iLpv v 2 Y V2 + e^u^, 
(6' 33 ) 2 = (2 + 5e l2 ^ c + 2e l4 ^ c )AV 2 ^z/ 2 

+ (1 + 4e i2t ^ c + e i4 ^ c )\v d e l ^u 3 Y U3 + e i2 ^ c v 2 d Y 2 3 , 
(b' 12 ) 2 = (2 + 5e^" c + 2e l4 ^" c )A 2 z/ 1 z/ 2 + (1 + e i2 ^ c + e i4 ^ c )Xv d e i(f ^ 2 Y ul 

+ ((1/2) + 2e i2 ^ c + (l/2)e i4 ^ c )A^e-^z/ 1 y; 2 + (1/2)(1 + e* 4 ^ 2 !^, 
(&' 13 ) 2 = (2 + 5e i2t ^ c + 2e i4 ^- c )A 2 z/iz/ 3 + (1 + e i2 ^ c + e^ c )\v d e llp " v z Y Vl 

+ ((1/2) + 2e i2 ^ c + (l/2)e i4 ^ c )A^e-^z/ 1 n 3 + (1/2)(1 + e i4 ^>X^ 3 , 
(6' 23 ) 2 = (2 + 5e^" c + 2e a ^ c )A 2 e i2 ^i/ 2 i/ 3 

+ ((1/2) + 2e i2 ^ + (l/2)e^)Au d e^(^n a + ^ 3 ) + e^^Y^ (B.5) 



e n = (4e i4 ^ c -4)A 2 z/ c ^e^ e - l ^z/ 1 y; i , 

e 22 = {2 - 2e l4 ^ c )X 2 u c v 3 u e i ^e i ^u 2 Y U2 , 

e 33 = (2 - 2e i4 ^)\W u e i ^e i ^iy 3 Y l/3 , 

e 12 = (2e i4 ^ c - 2)AV^e^V^z/ 2 Y; i + (1 - e^jAV^e^viy,,, 

e 13 = (2e* 4 ^ c - 2)A 2 z/ c ^e^e^z/ 3 ^ 1 + (1 - e^^X^vle^e'^v^, 

e 23 = (1 - e^AV^V^z^ + ^Y Vz ). (B.6) 



Let us discuss two particular limits where the formula becomes simple. In the limit 
M — > 00 and t> <j — > we obtain 



and 



(m eff )ij ~ Fi 




(B.7) 



where 




-2 



(B.8) 
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In the limit of vanishing phases, i.e. real VEVs, we obtain 



(m e ff\real) 



2 



(kv c + \v u v d )v° 




(B.9) 



where we have defined 



bi = Y Vi v d + 3Az/j. 



(B.10) 



Regarding the previous parameters we note that for the real case 




(B.11) 
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